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Abstract 

We say that a polygon inscribed in the circle is asymmetric if it contains no two antipodal 
points being the endpoints of a diameter. Given n diameters of a circle and a positive integer 
k < n, this paper addresses the problem of computing a maximum area asymmetric fc-gon 
having as vertices k < n endpoints of the given diameters. The study of this type of polygons is 
motivated by ethnomusiciological applications. 
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1 Introduction 


Imagine that we are at a concert of salsa and we want to retain the intrinsic rhythmic pattern 
to dance appropriately. Then we should know the clave Son. The clave Son rhythm might be 
represented as the 16-bit binary sequence 1001001000101000 or, as usual for cyclic rhythms, by 
onsets represented as points on a circle as in Figure [I] c). In this geometric setting, the clave Son 
is associated to a pentagon whose vertices are selected on a circular lattice of 16 points. So, the 
clave Son is represented by a special selection of k points on the circular lattice. This pattern 
has conquered our planet during the last half of the 20th century. But, what properties does this 
particular selection have? Notice that the polygon that represents the clave Son does not contain 
two antipodal points on the circle and moreover, it is easy to prove that this configuration is just 
the pentagon of maximum area without antipodal vertices (this later property produces a certain 
kind of asymmetry). Musicians have showed interest in finding similar patterns. Ethnomusicology 
is the discipline encompassing various approaches to the study of music that emphasize its cultural 
context. More specifically, Ethnomathematics is a domain consisting of the study of mathematical 
ideas shared by orally transmitted cultures. Such ideas are related to numeric, logic and spatial 
configurations mm. Related to spacial configurations, the area is a useful measure of evenness of 
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Figure 1: 
Cuba. 



a) Bulerfa b) Aka Pygmies rhythm c) Clave Son 

a) A flamenco rhythm in Spain, b) a rhythm used in Central Africa, c) the clave son in 



scales and rhythms in music theory [HIE]- Moreover, many musical traditions all over the world 
have asymmetric rhythmic patterns. For instance, the Aka Pygmies rhythmic pattern in Figure [l] 
b) has the so-called rhythmic oddity property discovered by Simha Arom [3j. A rhythm has the 
rhythmic oddity property if when represented on a circle it does not contain two onsets (the black 
points in the Figure [l]) that lie diametrically opposite to each other. Thus, the property asserts 
that one cannot break the circle into two parts of equal length, whatever the chosen breaking 
point, as there is always one unit lacking on one side. This property produces a kind of perceptual 
asymmetry. The asymmetry of the pattern is to some extent intrinsic, in the sense that there exists 
no breaking point giving two parts of equal length. Note that the oddity property requires that the 
circle is divided into an even number of units. The notion of rhythmic oddity has received different 
mathematical treatments; see m for more details. 

An algorithm for enumerating all the patterns satisfying the rhythmic oddity property has been 
proposed in j8]. Asymmetric rhythmics can also be found in the flamenco music of Spain (Figure 
[I] a)) and the clave Son in Cuba (Figure [I] c)). See [T0| [9j, and [16], for a detailed study on the 
preference of theses rhythms in their cultural contexts. 

Inspired in these ideas, we introduce the following geometric problem: 

Problem 1.1 Given n diameters in a circle and a positive integer k < n, select k endpoints of 
these diameters, no two of the same diameter, in such a way that their convex hull defines a k-gon 
of maximum area. 

Let us introduce some notation and related work. Let (po,Po), • • •, (p n -i,p' n -i) be n diameters of a 
given circle. Let S := {po,p'o,p\,p'\, ■ ■ ■ ,Pn-i,Pn-i} be an antipodal set containing the 2n endpoints 
of the given diameters. A sub-polygon (of S ) is a convex polygon whose vertex set is a subset of S. 
An antipodal polygon (on S) is a sub-polygon whose vertex set contains exactly one endpoint from 
each diameter (pi, p'fi of S [lj. An asymmetric polygon is a sub-polygon that contains no diameter. 
Therefore, an antipodal polygon is also an asymmetric polygon. Aichholzer et al. proved that an 
antipodal polygon of maximum area can be found in 0(n)-time [T]. The linear time algorithms 
they proposed are strongly based on a simple characterization for the extremal antipodal polygons. 
Namely, that an antipodal polygon of maximum area has an alternating configuration [1]. 

The problem however, is significantly different if we ask for an asymmetric A;-gon of maximum 
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area with k < n vertices in S. It is not difficult to come up with examples for which the simple 
characterization stated above does not hold if k < n. Aichholzer et al. presented an 0(n n ~ fc+1 )-time 
algorithm to compute an antipodal fc-gon of maximum area. However, the existence of a polynomial 
time algorithm to solve this problem was left as an open question [I]. 

In this paper, we answer this question in the affirmative. We distinguish two cases: If we are given 
a circular lattice with an antipodal set of 2n points (induced by n evenly spaced diameters), we 
show how to solve Problem |1.1| in constant time by providing a characterization of the solution. 
Otherwise, if the diameters are given in a general configuration, we show that the problem can be 
solved in 0(fcn 4 )-time using dynamic programming. 

The problem studied here is related to other optimization problems in mathematics. In compu¬ 
tational geometry, efficient algorithms have been proposed for computing extremal polygons with 
respect to several different properties [6]. Moreover, the so-called stabbing or transversal problems 
(see for instance [2j) belong to the same family that our problem. Recently, it has been proved 
that the following problem is NP-hard QH : Given a set S of line segments, compute the minimum 
or maximum area (perimeter) polygon P such that P stabs S, that is, at least one of the two 
endpoints of every segment s E S is contained in P. In operations research, global optimization 
techniques have been extensively studied to find convex polygons maximizing a given parameter [5]. 
Now, let us consider a different interpretation of the problem, let P be the convex hull of the given 
diameters. In this case, the solution to Problem 0 can be interpreted as the asymmetric polygon 
with k < n vertices that has its area closest to that of P. Therefore, this problem is related to the 
approximation of convex sets. In this setting, the best inner approximation of any convex set by a 
symmetric set is studied in m • Moreover, if we consider the “distance” to a symmetric set to be 
a measure of its symmetry [12], then our solution to Problem 0 provides the best approximation 
of a convex polygon inscribed in a circle by asymmetric sub-polygons. 

Finally, it is worth mentioning that although there exists a high number of applications of mathemat¬ 
ics to music theory, the research in music has illuminated problems that are appealing, nontrivial, 
and, in some cases, connected to deep mathematical questions. The problem introduced in this 
paper could be an example. 

The remainder of the paper is organized as follows. In Section [2] we consider the constrained version 
of Problem |1.1| in which the endpoints of the diameters are evenly spaced on the circle. In this 
case, we provide a characterization of the maximum k -gon which yields a constant time algorithm 
to solve Problem 0 The general version of this problem where the endpoints of the diameters are 


can be solved in polynomial time using dynamic programming. 


distributed anywhere on the circle is studied in Sectional In this case, we show that Problem 1.1 


2 Evenly spaced diameters 

In this section, we assume that S = {poiP'oiPhPi, ■ ■ ■ ,Pn-i,p' n -i} is the set of endpoints of n evenly 
spaced diameters on a unit circle. That is, S partition the circle into 2 n arcs of equal length. A 
£:-gon Qfc with vertices qo ,..., qk-i in S can be encoded by a sequence < ao,..., a*_i >, where a* 
is the number of arcs between the vertices qi and q%+\, 0 < i < k, qk = qo- For example, the clave 
Son in Figure [l] c) can be encoded by the sequence < 3,3,4, 2,4 >. In music theory, this sequence 
is called the interval vector (or full-interval vector) [14] . As with rhythmic patters, we assume that 
the sequence < ao, ■ ■ ■, a^-i > is cyclic, i.e., for i > k, = aj mo dfc. Throughout this paper, we 
identify a sub-polygon with k vertices by its corresponding encoding sequence. 
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Remark 2.1 Let Qbe a sub-polygon encoded by the sequence < ao,...,Ofc-i >■ The following 
properties hold: 


1. The area of the polygon Qk is - Jfi=o si n (^ a i 


2. The area is invariant under permutations of the sequence < ao, • • • ,Ofc-i >■ That is, for any 
0<i<j< k—1, the area of the polygons given by the sequences < ao,... ,a 2 ,... ,aj,..., ak- 1 > 
and < ao,..., a^,..., a*,..., a^-i > is the same. 


3. Polygon Qk is asymmetric if and only if there are no two indices i, j such that i ^ j and 
+ aj-|_i + • • • + aj = n. 


The following technical results will be used later. 

Lemma 2.2 Let a, b, n be positive integers with a — b > 2 and 2n > a + b. Then we have 

. ( 7ra\ . /vr6\ . (ir(a — 1)\ . / ir(b + 1) 

sin I — + sin — < sin - + sin - 

\ n J \ n J V n ) V n 


Proof. The previous inequality can be written as: 

2 sin ( cos ( *(*-<>) ) <2ain ( + cos / 2) 

V 2 n ) \ 2n ) \ 2n ) \ 2 n 

Because 0 < < 7r, 0 < 7r *- a ^~ 2 ^ < 7r ^~ b ^ < vr and from the fact that cos(-) is decreasing in 

(0,7r), this inequality holds. □ 


! /■ ^ ^ Lli 

For ease of notation, a sequence a,..., a of c elements all equal to a is denoted by The 

C 

following lemma characterizes the £:-gons (not necessarily asymmetric) of maximum area. 

Lemma 2.3 Given a regular polygon of m vertices inscribed in a circle, the maximum area sub¬ 
polygon with k < m vertices is encoded by any permutation of the following sequence 

< q+1, • • • ,q+ 1 , q, • • • ,q > , 

r k—r 


where m = kq + r, 0 < r < k. 


Proof. Let ao, ai,..., ak -1 be a sequence encoding a sub-polygon Qk with maximum area. Define 
^max — max (ao, ai,..., a^_r) and a m i n — min (ao, a \,..., a^—i). If a rnax Qmin 2, then we switch 


O'max to a ma 
O'max + ttmin 


- l. 


x mm 

and a m .i n to o. rn j n -{- 1. 


In this way, because k > 3 we have m = Jfi=o a * > 
and hence, Lemma 2.2 implies that the new polygon has larger area than the initial 
one leading to a contradiction. Therefore, a max ~ a. mm < 1. This implies that for some integer e 
and for every 0 < i < k — 1, either a* = e or a* = e + 1. 


Let u and v be the number of occurrences of the integers e and e + 1 respectively, that is, u(e) + 
v(e + 1) = m. Thus, we have {u + v)e + v = kq + r and ke + v = kq + r. Consequently, v — r is a 
multiple of k. Since v < k and r < k, we have only two cases: either v = r, or v = k and r = 0. In 
the former case, e = q and in the later, a* = e + 1 = q for every i. The result follows. □ 
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2.1 Case k odd 


In this section, we show that if the number of vertices k is odd, then there is a permutation of the 
sequence of Lemma 2.3 such that its corresponding sub-polygon is asymmetric. 


Theorem 2.4 Given a regular polygon of 2 n vertices inscribed in a circle and positive integers 
k, q, r with 2n = kq + r,0<r< k, the maximum area asymmetric sub-polygon with k vertices is 
given by the following sequences: 

For r even, 

<<7+1, ••• ,<7+1, < 7 , ••• ,< 7 ,g+l, ••• ,<7 + 1, <7, ••• ,q> 

'-v-•' '- v -•' ^- v -^- v -•' 

r/2 (k—r— 1)/2 r/2 (k—r+l)/2 

For r odd, 


< q, 


,<7, Q+ 1, 


,q + 1, q, 


,q, <7 + i, 


> q +1 > 


(k—r)/2 


(r- 1)/2 


(k—r)/2 


(r+l)/2 


Proof. We note first that both sequences are composed by r occurrences of the element <7 + 1 and 


k — r occurrences of the element q. Therefore, both fc-gons are maximum area by Lemma 2.3 
-i) for 0 < i < , we say that a. + (fc-i) is the homologous element of a* in 


Moreover, a* = a (& 


the sequence and vice versa. Notice that, except for aj--i, every element in the sequence has its 
homologous element. 


We now prove that the encoded polygon is asymmetric. Let A and B be a partition into two 
subsequences of the proposed cyclic sequence so that A contains consecutive elements and B 
the other elements. 

Let Sa and Sb be the sum of the elements of A and B , respectively. If a^_i ^ A, then B consists 
of the — 9 ^ homologous elements of the corresponding of A and a^-i and hence Sb = Sa + a-k-i- 
On the other hand, if ak-i G A, then B contains homologous elements of the ones in A plus 

two more elements, say ai,aj (recall that ak-i has no homologous element). Therefore Sb — Sa = 
ai + aj — ak-i for some indices i j. Regardless of the case, we have that Sb > Sa, Sa < n and 
Sb > n. As a consequence, any sum of 2 ^ consecutive elements is different from n. 

For an arbitrary partition of the sequence, we proceed in a similar way. If the sizes of A and B are 
arbitrary, we also have two cases. If |A| < then Sa < n and if |L>| < then Sb < n. In 

both cases, the sequences correspond to asymmetric polygons. 


□ 


2.2 Case k even 


Unfortunately, if the number of vertices k is even, the following result prevents us from seeking the 
solution as a permutation of the sequence of Lemma |2.3[ 

Theorem 2.5 If k is even, then no permutation of the sequence 

< <7 + i, • • • ,q + 1, <7, • • • ,q> 

<_ v _< S. - ^ 

r k-r 


corresponds to an asymmetric polygon. 
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Proof. Set f(i) = a-i + aj+i + ... + a i+k / 2 _ i- We have f(i + 1) — f(i ) = a i+k / 2 — a*, which implies 
that f(i + 1) — /(*) E {—1,0,1}. If /(0) = n, then the sequence is non-asymmetric. Assume that 
/(0) < n. Thus, f(k/2) > n (/(0) + /(fc/2) = 2n). Consider now the discrete function given by 
the integer values /(0),/(l),..., /(fc/2). Since \ f(i + 1) — f(i)\ is equal to 1 or 0, /(0) < n and 
/(fc/2) > n, it follows that there exists an index j, 0 < j < fc/2, such that f(j) = n and the 
sequence generates a non-asymmetric polygon. The proof for the case /(0) > n is analogous. □ 

In the following, we prove that it is possible to find the solution for the case fc even by considering 
the sequence that provides the second maximum area. 

Theorem 2.6 If k is even, then any permutation of the sequence 

<q + l, • • • ,q+ 1 , q- 1 , g, • • • ,q> 

S. - ^ ^ ^ v _✓ 

r+1 fc— r— 2 

corresponds to the sub-polygon of k-vertices with the second maximum area. 


Proof. We first prove that if a* is an element of a sequence providing the fc-gon with the second 
maximum area then a* E {q — 1, q, q + 1, q + 2} for i = 0, ..., fc — 1. 

Denote by a max and a m i n the values maxjao, a ±,..., afc-i} and minjao, ai,..., afc-i}, respectively. 
Suppose that a max > q + 3. Note that if a m i n > q + 1, then X)i=o a * > which is a contradiction. 
Therefore, we assume that a m i n < q. In this case, switching a max by a max — 1 > q + 2 and a m i n by 
o-min + 1 yields a larger polygon according to Lemma [272] However, the new polygon cannot be the 
largest one because it contains an element aj > q + 2, in contradiction with Lemma 

Suppose that a m i n < q — 2. Following a similar argument to the above case, a max > q + 1 and by 
switching a max and a m * n by a max — 1 and a m i n + 1 < q — 1, respectively, the area of the new polygon 
increases. It cannot be the largest one because it contains an element aj < q — 1. 

Using similar arguments we can show that the elements <7 + 2 and q — 1 appear at most once in the 
fc-gon with the second maximum area. 

Let A, B , C, D be the number of occurrences of q — 1, q, q + 1, q + 2, respectively. The we have: 


2.3 


A+B+C+D 
A(q — 1) + B(q) + C(c[ + 1) + D[q + 2) 
(A + B + C + D)q + {-A + C + 2D) 

-A + C + 2D 
C 


k 

2 n 

kq + r 
r 

r + A- 2D 


Now, the area of the fc-gon we are looking for is: 


S(A, B, C, D) = A sin 


7T(<7 - 1) 


n 


+ B sin 


irq 

n 


+ C sin 


Tr(q + 1) 


n 


+ D sin 


tt(<7 + 2) 


n 


Since 0 < A,D < 1, we have four cases to consider. Case (A, D) = (0,0) corresponds with the 
fc-gon of maximum area (Lemma 2.3). Hence, the areas of the remaining cases are: 


5(1, Hi, Ci, 1) = sin (^ 7r(g ~ 1) ) + (fc - r - 1) sin (^) + (r - 1) sin () + sin (^ 7r(9 ^~ 2) ) 

5(0, f? 2 , C 2 ,1) = (fc - r + 1) sin (^) + (r - 2) sin + sin ( 7r(g + 2) ) 

5(1, B 3 , C 3 , 0) = sin ( 7r(g ~ 1) ) + (fc - r - 2) sin (^) + (r + 1) sin j 
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We claim that 5(1, B\, C\, 1) < 5(0,52,(72,1) < 5(1,53,(73,0). To prove this claim, note that 


5(0, 5 2 , (7 2 ,1) - 5(1, 5 l5 (7i, 1) = 2 sin ™ 


vr(g - 1) 


77- 


— sm 


?r(g+ 1) 


n 


By using Lemma 2.2 with a = 
Notice that 2 n > kq > 4q > 2q 


q + l,b = q — lwe have that 5(0, 5 2 , (7 2 ,1) — 5(1, 5i, (7i, 1) > 0. 
= a + b and a — b = 2 in this case. 


To prove the other inequality, we use the following trigonometric identities: 


sin(a) — sin(/3) = 2cos((a + /?)/2) sin((a — /3)/2) 
sin(3a) = 3sin(a) — 4 sin 3 (a) 

Using these identities, we obtain the following: 


5(1,5 3 , (7 3 ,0) — 5(0,5 2 , (7 2 ,1) = sin 

= 6 cos 

= 2 cos 

= 2 cos 


n ^- 1 ^ + 3 a uJ 7T(q + l1 


n 

ir(2q + 1) 
2 n 

vr(2 q + 1) 
2 n 

7r(2g + 1) 
2 n 


n 


7r 


sin ( — ) — 2 cos 
2 n 


n q 
n 

n(2q + 1) 


— 3 sin ( — ) — sin 
n 


ir(q + 2) 


n 


2 n 


sm 


37t\ 

2 n) 


3 sin 


4 sin d — 
2 n 


7T 

2n 

7T 


— sm 


3-7T 
2 n 


Because 2?r = /eg + r, > 4 and from the fact that (7 2 = r — 2 > 0, we get that 2n > 4q + 2 and 
f ■ As a consequence, 5(1, 5 3 , (7 3 ,0) — 5(0, 5 2 , (7 2 ,1) > 0 yielding the result. □ 


The above result asserts that any sequence with coefficients A = 1, B = k — r — 2, C = r + 1, D = 0 
provides the fc-gon with the second maximum area. In the following result, we provide a sequence 
that in addition corresponds to an asymmetric £:-gon. 

Theorem 2.7 The sequence corresponding to the maximum area asymmetric k-gon for k even is 
as follows: 


<q - i, q + i, 


,q + 1, 


,q, q + i, 


,q + i, q, 


,q> 


r/2 


( fe — r — 2)/2 


(r+2)/2 


(k—r—2)/2 


Proof. The proof is similar to that of Theorem |2.4| The sequence proposed, say Q : is a permutation 
of the sequence in Theorem 2.6 Therefore, its corresponding k -gon has the second maximum area. 
We show that it is also asymmetric. 


Observe that for 0 < i < /c/2, it holds that a* = a i+k / 2 ■ We say that a i+k / 2 and a* are homologous 
elements for 0 < i < k/2. 

Let A and 5 be a partition into two subsequences of the cyclic sequence Q. First, suppose that A 
(and 5) consists of k/2 consecutive elements. In this case, neither A nor 5 contains homologous 
elements. Let Sa and Sb be the sum of the elements of A and 5, respectively. If ao £ 4 then 
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a k /2 £ B and Sa~ Sb = clq — a k / 2 = —2 / 0. Similarly, if ao £ B we have Sa — Sb / 0 proving 
onr claim. 

Suppose now that the subsequence A contains less than k/2 consecutive elements. Then, every 
element in A has its homologous element in B except maybe for oo or a k / 2 ■ If ao £ A, then 
Sb ~ Sa = a k /2 — ao + a u + a v + ..., where a u ,a v ,... are elements of B having no homologous 
element in A. In this case, Sb > Sa- The case \B\ < k/2 is analogous. Hence, there are no indices 
i,j, i / j. so that a* + aj+i + ... + a,j = n. This completes the proof of the theorem. □ 


3 The general case 

Dfaz-Banez et al. m showed that the following problem is proved to be NP-hard: Given a set 
of n line segments in the plane (or in a circle), find a maximum-area convex polygon having as 
vertices at least one endpoint from each segment. They also show that if the segments are pairwise 
disjoint, then the problem can be solved in polynomial time. A constrained version in which each 
segment joins two antipodal points on a circle was recently studied by Aicholzer et al. pQ. They 
showed that this constrained problem can be solved with a simple linear time algorithm and they 
asked about the hardness of the problem of selecting k endpoints instead of n. As a consequence 
of the hardness result of Diaz-Banez et al. m, this problem becomes NP-hard if we remove the 
antipodality constraint. 

A dynamic programming algorithm presented in this section finds the asymmetric fc-gon of maximum 
area in time 0(fcn 4 ). For the special cases k = 3 and k = 4 we show how to solve the problem in 
linear time. 


3.1 Dynamic programming algorithm 


Let S := {po,PoiPhPii ■ ■ ■ )Pn-i,Pn_i} be the set of the endpoints of n diameters in the unit circle 
and let O be the center of this circle. We define a wedge ( i,j ) as a convex polygon being the convex 
hull of O and a subsequence of points (need not be consecutive) of S visited when walking clockwise 
along the circle from pi to pj, see Figure 2a We say that pi and pj are anchors of a wedge (i, j). 




Pi 



Figure 2: The figure a) shows a wedge with anchors pt and pj, b) and c) show the two possible 
configurations of a double-wedge ( i,j,l,t ). 


A double-wedge ( i,j,l,t ) is the union of two disjoint wedges (i,j) and ( l,t ) such that: (1) this union 







Figure 3: A maximum area asymmetric k-gon on S formed by the union of a double-wedge ( i,j , l, t) 
and the triangles Opjpi and Opipt . 


contains no diameter and (2) either p' t lies clockwise between pi and pj , or p'j lies clockwise between 
pi and pt ; see Figure [2] for an illustration. The anchors of a double-wedge are the anchors of the 
wedges that form it. 

Lemma 3.1 The maximum area asymmetric k-gon (k > 3) on S is formed by the union of a 
double-wedge ( i,j,l,t ) and the triangles Opjpi and Opipt- 

Proof. Let Q be the maximum area asymmetric k- gon on S. Let pi be a vertex of Q and let p[ be 
its antipodal point. Let pj ( resp. pi) be the first vertex of Q found when walking counterclockwise 
(resp. clockwise) from p[ along the boundary of the circle. Let pt be the first vertex of Q found 
when walking counterclockwise from pi along the boundary of the circle. See Figure [3} 

We prove by contradiction that pj , pi and pt exist and are different from pi. Let Sq C S be the 
set of vertices of Q and its antipodal points and the asymmetric fc-gon Q of maximum area on S 
is a maximum area antipodal on Sq. Suppose that pj does not exist or matches with pi, then all 
vertices of Q lie in the same semicircle of the two defined by the diameter pip\, then Q would not 
be a maximum area antipodal polygon on SqJT , this is a contradiction. The proof that pi and p t 
exist and are different from pi is analogous, note that pi can matches with pt. The union of the 
wedges (i, j) and ( l,t ) contains no diameter because Q is an asymmetric polygon. We now prove 
that p'j lies clockwise between pi and pt or p' t lies clockwise between pi and pj. If p'j lies clockwise 
between p t and pi then p' t lies clockwise between pi and pj, if p'j lies clockwise between pi and pt we 
are done, finally if p'j lies clockwise between p[ and pi we obtain a contradiction because all vertices 
of Q would lie in the same semicircle of the two defined by the diameter pjp'j. Therefore, Q is not 
a maximum area antipodal polygon on Sq [I] . □ 


When a double-wedge uses s points of S we say it is an s-double-wedge. For 3 < s < k, let 
f(i,j , l , t, s ) be the area of the s-double-wedge of maximum area with anchors (■ i,j , l, t ). Notice that 


for s = 3 either pi and pj , or pi and pt coincide. By Lemma 3.1, the area of the maximum area 
asymmetric k- gon is equal to the maximum of {f(i, j, l,t,k ) + Ao PjPl + Ao PiPt }, where (i. j. I, t) goes 
over all the possible anchors of a /c-double-wedge. Therefore, having computed f(i,j,l,t,k ) for all 
possible fc-double-wedges, the maximum area asymmetric fc-gon can be found in time 0(n 4 ). 


Lemma 3.2 All possible values of f(i,j,l,t,s) with 3 < s < k can be computed in time 0(kn 5 
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Proof. For s = 3. if i = j then f(i, j, l,t,3) = Ao PlPt else f(i, j, l,t,3) = Ao PiPj ■ All possible values 
of 3) can be computed in time 0(?i 3 ). For fixed values i, j, l, t and s (s > 3), the value 

of f(i,j , l, t, s ) can be computed in time 0(n). If p' t lies clockwise between pt and pj then: 

f(i,j,l,t,s ) = ma x{f(i,m,l,t,s - 1) + A 0pmPj } 

for each p m between pi and pj in clockwise, p m ^ p[ and p m / p[. Otherwise, if p'j lies clockwise 
between pi and pt then: 


f(i,j,l,t,s) = ma x{f{i,j,l,m,s - 1) + A 0pmPt } 

for each p m between pi and pt in clockwise, p m / p[ and p m / p'j. In fact, all possible values 
of ) with a fixed value s (s > 3) can be computed in time 0(n 5 ) having computed the 

values f(i,j,l,t,s — 1). Then, by iterating s from 3 to k, all possible values of f(i,j,l,t,s) can be 
computed in time 0(kn 5 ). □ 


We now show how to improve this running time by computing all possible values of l, t, s ) for 

a fixed value s > 3 in time 0(?r 4 ). The following lemma is depicted in Figure 4a 



Lemma 3.3 Let f(i,ji,l,t,s ) and f(i,j 2 ,l,t,s) be the areas of the maximum area s-double-wedges 
with anchors (i,ji,l,t) and ( i,j 2 ,l,t ), respectively, such that s > 3 and p\ lies clockwise between 
Pi and pj x . Let Pi, ■ ■ ■ ,p r ,Pj i and Pi, ■ ■ ■ ,p q ,Pj 2 be the subsequences in S of the wedges (i,ji) and 
(i, j' 2 ), respectively. If pj 1 precedes pj 2 when walking clockwise along the circle from pi, then p r also 
precedes p q . 

Proof. If p q lies clockwise between pj 1 (including it) and pj 2 we are done. Now, consider the case 
when p q lies clockwise between pi and pj 1 . The subsequences Pi, ■ ■ ■ ,p r ^Pj! and pt, ■ ■ ■ ,p q ,Pj 2 are 
part of the maximum area s-double-wedges (i,ji, l, t) and [i,j 2 , l, t), respectively. Therefore we have 

f{i,ji,l,t,s) = f(i,r,l, t,s — 1 ) + A 0prPjl 

> f(i , m, l, t, s — 1) + Ao PmPn for each p m lies clockwise between pt and p 3] (1) 
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for each p rn lies clockwise between pi and pj 2 (2) 


f(i,j 2 ,l,t,s) = f(i,q,l,t,s- 1 ) + A 0pqPj2 

> f('i, rn, l, t, s - 1) + A 0prnPj2 

From [l] and [2] we obtain the following inequalities: 

f(i, r, l, t, s - 1) + A 0prPjl > f(i, q, l, t, s - 1) + A 0pqPjl 

f(i, r, l, t, s - 1) + A 0prPj2 < f(i, q, l,t,s- 1) + A 0pqPj2 
adding both we obtain that 

Ao PrPj2 ~ Ao PrPjl < Ao PqPj2 — Ao PqPjl (3) 


Suppose that p q and p r are inverted. Let z, x and y be the angles between p q , p r , pj t and pj 2 (see 
Figure [5]), then Inequality [3] can be written as: 

sin(x + y) — sin 2 ; < sin(a: + y + z) — sin(^ + x) 


o ( ( i_ , y\ • (V\ 

2 cos(x + -) sin(-) < 2 cos(x + z + -) sin(-) 


/ - y \ ^ ( - - 2 / \ 

cos(x + -) < cos(x + z + -) 


( 4 ) 


Because z+x+y < it and cos(-) is decreasing in the interval (0, 7 r), Inequality [4] yields a contradiction. 



Figure 5: This figure shows the triangles involved in the proof of Theorem 


3.3 


□ 


The following lemma has an analogous proof to that of Lemma 3.3 and is depicted in Figure [4b 


Lemma 3.4 Let ) and f(i,j,l,t 2 ,s ) be the areas of the maximum area s-double-wedges 

with anchors (i,j,l,t\) and (i,j,l,t 2 ), respectively, such that s > 3 and p'j lies clockwise between 
Pi and pt x . Let pi,... ,p r ,pt q and pi,... ,p q ,pt 2 be the subsequences in S of the wedges ( l,t \) and 
(l,t 2 ), respectively. If p^ precedes pt 2 when walking clockwise along the circle from pi, then p r also 
precedes p q . 


Theorem 3.5 Let S := {poiPoiPhPii ■ ■ ■ >Pn-i,p' n -i} be the set of the endpoints of n diameters in 
the unit circle and let k be an integer such that k < n. The asymmetric k-gon of maximum area 
can be found in time 0(kn 4 ). 
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Proof. Having computed l, t, k) for all possible /c-double-wedge, the maximum area asymmet¬ 
ric fc-gon can be found in time 0(n 4 ). We show how to compute f(i,j , l, t, s) for all s-double-wedge 
with 3 < s < k in time 0[kn 4 ). 

For s = 3, 3) can be computed in time 0(n 3 ). For s > 3, if p' t lies clockwise between pi 

and pj , f(i, j, l, t, s ) = max{/(«, m, l,t,s — \) + A Op mPj }, for each p m between pi and pj in clockwise, 
Pm Pi and p m ^ p\. Suppose that we have computed for i, j, l and t the value r such that 
, t, s ) = f(i, r,l,t,s — 1) +Ao PrPj ■ Then, to compute f(i,j + 1, l, t, s ) we do not need test each 
p m between pi and Pj+i, we need only to test each p m in the interval from p r to Pj+i ■ In fact, we 
can compute ) for all pj that follows p' t in clockwise for fixed values i, l, t and s in time 

0(n). The other case, when jj'- lies clockwise between pi and pt is analogous. In fact, all possible 
values of l , t, s ) with a hxed value s (s > 3) can be computed in time 0(n 4 ) having computed 
the values f(i,j,l,t,s — 1). Then, iterating s from 3 to k, all possible values of f(i,j,l,t,s ) can be 
computed in time 0(kn 4 ). □ 

3.2 Special case k = 3 

Let S := {po,p' 0 ,pi,p'i,... ,Pn-i,p' n -i} be the set of endpoints of n diameters in the unit circle, 
7i > 3. In the following, we prove that the triangle of maximum area on S contains no diameter. 
Using this characterization, we show how to find the largest triangle in linear time. We use pipj to 
denote the clockwise arc from pi to pj and Zpipjpf. to denote the inner angle of the triangle PiPjPk 
with apex p 3 . 

Lemma 3.6 A maximal triangle with vertices in S contains no diameter. 

Proof. Suppose for a contradiction that the maximal triangle contains a diameter. Let PiPjp'i be 
the maximal triangle. Let pi and p\ be another diameter of points in S. If pi is in p t pj then p\ is 

in p'ip'j. The triangles PiPjp'i and PiPjp\ have a common base, the segment Pip'j, the height of PiPjp'i 
from p[ is smaller than the height of PiPjp\ from p\\ see Figure 0 Therefore, A pip . p , > A pjP;/p , which 

is a contradiction. The proof in the case pi is in pjp\ is analogous. □ 



Figure 6: A p ^ pjPi > A p . p i p . because the height of p\pjpi from p[ is greater than the height of Pip'iPj 
from pi. 

Note that each triangle contains either two angles greater or equal than 2-/r/3, or two angles smaller 
or equal than 27r/3; we call such pair of angles the base angles of the triangle. A critical vertex is 
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a vertex supporting a base angle. The following lemma is key to our algorithm. 

Lemma 3.7 Let c\, C 2 , p, 62 and b\ be points in the unit circle that lie in this order when walking 
counterclockwise along the circle from c\. 


1. If Zpb\Ci < 27t/3 and Apc\b\ < 2ir/3, then the maximal triangle that contains p is c\pb\ 

2. If Zpb 2 C 2 > 27t/3 and Zpc 2&2 > 27r/3, then the maximal triangle that contains p is C 2 pb 2 


Proof. Using similar triangles we can show the following inequalities for the first case (see Fig¬ 
ure 7a): 

Ap C2 b‘2 A -Apc 2 bl A Ap C 1 b I 
h 2 ^ Ape 1 1)2 A Apr i }) | , 


and these ones for the other case (see Figure [7b]): 


A, 


pC2b2 


> A 


pc2bi 


> A, 


pcibi 


A 


pc 2 b 2 


> A. 


pcib 2 


> A 


■pcibi 


□ 


V 



Figure 7: This figure shows triangles with p as 
is the closest to the equilateral triangle pbc. 



■itical vertex. In both cases the maximal triangle 


Theorem 3.8 Let S := {poiP'oiPhPii ■ ■ ■ ,Pn-i,Pn-i} be the set of endpoints of n diameters in the 
unit circle, n > 3. The triangle of maximum area with vertices in S can be found in linear time. 

Proof. The idea of our algorithms is to iterate among the endpoints of the n diameters in clockwise 
order computing for each pi the maximal triangle with pi as critical vertex. Let M be the triangle 
of maximum area with vertices in S. Then M = max{T,;}, where Tj is the maximal triangle with pi 


13 








Pi 



Figure 8: The two candidates to be the maximal triangles with critical vertex pi, PiP^pPm, and 
PiPkiPli- 


as critical vertex. Lemma 3.7 allows us to find T % comparing just two triangles PiPj,Prm and PiPkiPu 
(see Figure [8]), where j t < ki < U < rrii and the following inequalities hold: 


PiPj, < 2vr/3 <KPj(+i (5) 

PiPki > 2tt/3 >Pi'Pk-- \ (6) 

pipii< 4tt/3 <ppp^ x (7) 

Piihm > 4vt/3 >ppp^Z\ (8) 


Note that ji + 1 > ji, k l+ i > ki, U + \ > k and rri l+ \ > rrii. If we have for a point pi the points p 3i , 
Pki, Pli and Prm, then moving from these points in clockwise we can find for pi+i the points Pj i+1 , 
Pk i+ 1 , Pi i+1 and Pm i+ 1 - Therefore, we can compute the maximal triangle with p t as critical vertex 
for all pi in S in linear time. □ 


3.3 Special case k — 4 

Unlike the triangle, the quadrilateral of maximum area could contain diameters. However, we prove 
that the asymmetric quadrilateral of maximum area satisfies a property that allows us to find it in 
linear time. 

Lemma 3.9 Let S := {po,PoiPi,Pi, ■ ■ ■ ,Pn-i,p' n -i} be the set of endpoints of n diameters in the 
unit circle, n > 4. Let pi, pi, p’, and p[ be the vertices of a maximal quadrilateral. Then the arcs 

p#j, PlPt, rtp'j and ifp' t contain no other points of S, see Fz 5 we[i] 

Proof. Certainly pi is the farthest point of S from Pi;p' } in pi;p' } or p' t is the farthest point of S from 

Pip'j in p'jPi- Suppose that pi ( p' t ) is the farthest point of S from pip'j in pppj (p'jPi), then p[ (p t ) is 

the farthest point of S from pip'j in pLp t (pip'j) . In fact p' t ( pi ) is the second farthest point of S from 

Pip'j in p'jPi ( pip'j ). Hence the points pi and pt (p[ and p{) are consecutive in clockwise. Analogously 
we can prove that pi and pj {p[ and pb) are consecutive . □ 
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Note that the four possible asymmetric quadrilateral with vertices in the subset { pt , pj , pi, pt, p[, 
p'j, p'j. p' t } have the same area because the segments pipb and p[pj ( pip[ and p[pt) are parallel and 
equal. As a consequence of Lemma 3.9 pj = Pi+i and pt = pi+i- From the above lemma we infer 
the following corollary. 


Corollary 3.10 Let S := {po,p'o,Pi,p'i, ■ ■ ■ ,p n -i,p' n -i} be the set of endpoints of n diameters in 
the unit circle, n > 4. Let Q be a maximal asymmetric quadrilateral on S. If pi £ S is a vertex of 
Q then its opposite vertex in Q is p' i+l or p , i _ 1 . 

Lemma 3.11 Let a, b, c, d, a', e, / and d! be clockwise points in a circle, where aa' and dd! are 
diameters. If d is farther than c from ae, then d is farther than c from bf. 


Proof. Let m be the intersection point of bd and cf, let n be the intersection point of ad and ce. 
Because d is farther than c from ae, A ace < A a d e . Because triangle ane is common for the triangles 
ace and ade , A anc < Ad ne - Since the triangles anc and dne are similar, ac < de. Moreover, be < ac 
and de < df as they are in the same semicircle, respectively. From these inequalities, we infer that 
be < df. Because the triangles bmc and dmf are similar, we know that Ab mc < Admf- Since triangle 
bmf is common for the triangles bef and bdf, A^ c f < Therefore, d is farther than c from bf ; 

see Figure [lOj □ 

We are now ready to show the main result of this section. 

Theorem 3.12 Let S := {poiPoiPhPii ■ ■ ■ ,Pn-i,Pn-v} be the set of endpoints of n diameters in 
the unit circle, n > 4. The asymmetric quadrilateral of maximum area can be found in linear time. 


Proof. Let Q be the asymmetric quadrilateral of maximum area on S and let Q, = PiPip' i+ iP'i +1 
a maximal asymmetric quadrilateral with Pip' i+ i as a diagonal, using Corollary 3.10 we have that 
Q = max{(5i}. The idea of our algorithm is to iterate among the endpoints of the n diameters in 
clockwise order and compute, for each pi, the maximal asymmetric quadrilateral Qi with Pip [ + 1 as a 
diagonal. Suppose that for an arbitrary point pi E S we have computed Qi = PiPip'i + iP'i + i w here Pi 
and p\ + i are the farthest points (the farthest and the second farthest maintaining the asymmetry) 
from Pip' i+ i on S. To compute Qi+i we need to find p m such that p m and p' m+ i are the farthest 
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Figure 10: Representation of lemma 


3.11 


we can infer that p m and p' m+ \ lie on p/p' +2 and 

p[ +1 Pi- )-i, respectively (see Figure |ll|. Therefore, making a clockwise iteration on S we can find 
all the maximal asymmetric quadrilaterals in linear time, computing for each pj the farthest points 
(the farthest and the second farthest maintaining the asymmetry) pi and pi+\ from Pip[ + i on S 
using the farthest points found in the previous step of the iteration. □ 


points from pi + ip' i+2 on S , using Lemma 




3.11 



Figure 11: pi is the farthest point of S from p*p( +1 and p rn is the farthest point of S from pi + ip' i+2 . 
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